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Abstract 

Autism Spectrum Disorders (ASD) are complex neurodevelopmental diseases clinically defined by dysfunction of social 
interaction. Dysregulation of cellular calcium homeostasis might be involved in ASD pathogenesis, and genes coding for the 
L-type calcium channel subunits Ca v 1.2 {CACNA1C) and Ca v B2 {CACNB2) were recently identified as risk loci for psychiatric 
diseases. Here, we present three rare missense mutations of CACNB2 (G167S, S197F, and F240L) found in ASD-affected 
families, two of them described here for the first time (G167S and F240L). All these mutations affect highly conserved 
regions while being absent in a sample of ethnically matched controls. We suggest the mutations to be of physiological 
relevance since they modulate whole-cell Ba 2+ currents through calcium channels when expressed in a recombinant system 
(HEK-293 cells). Two mutations displayed significantly decelerated time-dependent inactivation as well as increased 
sensitivity of voltage-dependent inactivation. In contrast, the third mutation (F240L) showed significantly accelerated time- 
dependent inactivation. By altering the kinetic parameters, the mutations are reminiscent of the CACNA1C mutation causing 
Timothy Syndrome, a Mendelian disease presenting with ASD. In conclusion, the results of our first-time biophysical 
characterization of these three rare CACNB2 missense mutations identified in ASD patients support the hypothesis that 
calcium channel dysfunction may contribute to autism. 
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Introduction 

Autism spectrum disorder (ASD) is defined by dysfunction of 
social interaction and communication, stereotypic behavior and 
sensory integration problems. It is a complex neuro-developmental 
disease, which might result from an altered brain ontogenesis or 
altered neural homeostasis. According to the Centers for Disease 
Control and Prevention (CDC) the prevalence of ASD in the US 
population is 1:88 [1]. Until now, genetic explanations for ASD 
are limited to rare chromosomal abnormalities like copy number 
variations or very rare single gene disorders. A first hint to unravel 
pathophysiological pathway of ASD came from the identification 
of the mutation p.G406R found in the L-type calcium channel 
pore-forming subunit (Ca v 1.2) gene CACMA1C in patients with 
Timothy Syndrome (TS) [2]. Whole-cell patch clamp recordings 
showed that the TS-mutation leads to a decelerated and 
incomplete inactivation of the calcium inward current. The 
phenotype of TS demonstrates the consequences of inadequate 
inactivation behavior of voltage gated calcium channels (VGCC) 
in different biological contexts, like heart, brain, and the immune 
system. Calcium channel inactivation is a key mechanism by 
which cells are able to tighdy control intracellular calcium levels 



and therefore the activity of excitable cells. Of note, calcium 
channel inactivation is not exclusively controlled by the channel 
pore, but also depends on auxiliary subunits, namely Ca v R and 
Ca v ot28 [3,4]. The assembly of the subunits and their specific 
isoforms determines the distinctive behavior of VGCCs for the 
neuronal function. The number of permutations of calcium 
channel complexes with their subunits, isoforms and splice- 
variants paves the way for fine-tuned calcium channel function 
adapted for information processing. There are transcripts of ten 
pore-forming Ca v Q(l -subunit genes (CACNA1 A to CACNA1I and 
CACNA1S), which undergo extensive alternative splicing and in 
concert with auxiliary subunits exhibit different biophysical 
properties and expression profiles as well as distinct subcellular 
targeting [5]. 

Particularly Cav(3-subunits have been shown to impact surface 
expression and modulation of channel activity and kinetics, 
leading to an increased L-type calcium channel activity, as 
revealed in whole-cell [6,7] and single-channel recordings [8]. A 
wide variation of inactivation behavior has been described for the 
splice variants of the B2-subunit [4]. Even slight structural 
differences within the CavP-subunit can strongly modify the 
gating behavior of L-type calcium channels [8-10]. Interestingly, 
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the expression of distinct CavP-subunits within the brain is 
dependent on the stage of neuronal development [1 1] . Since the 
discovery of the TS-mutation and its effect on inactivation of the 
L-type calcium channel, additional calcium channel genes or loci 
have been associated with ASD, e.g. the genes for the pore- 
forming subunits CACNA1D [12], CACNA1F [13], CACNA1G [14], 
CACNA1H [15] and the auxiliary subunit CACNB1 [14]. These 
studies emphasize the importance of calcium channel activity, 
inactivation and calcium signaling in a broader sense for the 
pathophysiology of ASD. Interestingly, the activity and kinetics of 
several pore-forming subunits of L-type (Ca v l. 1-1.4) and Non-L- 
type (Ca v 2. 1-2.3) channels [2,12,13,16] are regulated by CavP- 
subunits. Moreover, they are involved in various signaling 
pathways [2,17,18], all of which have been linked to ASD 
pathophysiology previously. Nevertheless, studies considering the 
role of calcium channels in ASD have focused on the pore-forming 
subunits of the calcium channel complex, even though gene 
clusters of interacting proteins participating in linear signaling 
pathways would have a similar chance of being involved in the 
etiology of ASD [19]. 

Based on our previous studies showing that current and gating 
kinetics profoundly depend on the particular Ca v P-subunit 
isoform or splice variant associated with the channel pore [8,9], 
we hypothesized that Ca v P-subunit mutations might lead to 
electrophysiological phenotypes similar to that observed in TS. 
Though all four Ca v P-subunit genes (CACNB1-CACNB4) are 
expressed within the brain [20] we chose CACNB2 because of 
positional evidence from a meta-analysis of linkage data [21]. 
Here, Trikalinos and colleagues showed genome-wide suggestive 
significance for a designated bin 1 Op 1 2-q 11.1, which embraces a 
large genomic region including the CACMB2 gene in autistic 
siblings. Recently, CACMB2 was found as a risk locus for five major 
psychiatric disorders including ASD [22] and thus is regarded as a 
susceptibility gene. 

In our current study, we searched for mutations in the P2- 
subunit gene that might affect gating behavior of voltage- 
dependent calcium channels. In a mutation screening, 155 patients 
with ASD were sequenced and the results were filtered for most 
promising candidates. We included only missense variants, which 
were not present in our cohort of 259 matched controls and 
located in highly conserved regions of the protein as an indicator 
of functional importance. Furthermore, the mutations had to be 
unknown or potentially damaging according to in-silico predic- 
tions. 

Here, for the first time, we present three missense mutations 
located in conserved regions of the calcium channel CACMB2 gene 
found in three families affected by ASD. In electrophysiological 



analyses of recombinant channels, these missense mutations were 
found to differentially alter current kinetics. Though our data do 
not prove an association between these CACNB2 mutations and 
ASD, our findings support the idea of CavP2 variants being of 
functional relevance for ASD pathophysiology. 

Materials and Methods 

Ethics Statement 

Procedures were approved by the Institutional Review Board 
(application number 04-223) of the medical faculty of the 
University of Cologne. 

Subjects 

gDNA samples were available from 259 healthy controls and 
155 patients with Autistic Spectrum Disorder (ASD), each are 
ethnically matched groups of Caucasians of west-European origin. 
The cases were elicited from three different sources: Twenty 
patients were recruited by the Department of Child & Adolescent 
Psychiatry and Psychotherapy, University of Cologne. These 
patients meet diagnostic criteria for 'broad ASD' (e.g. Aspergers's) 
based on interview for DSM-IV-R. A written informed consent 
was obtained from the parents or legal guardians of the minors 
using the consent form approved by the I.R.B. of the University of 
Cologne (Germany). Ten cases were included from Autism Tissue 
Program (www. atpportal.org) from the Maryland National 
Institute of Child Health and Human Development (NICHD) 
Brain Tissue Center and the Harvard Brain Tissue Resource 
Center (USA), published e.g. in [23]. 125 cases were provided by 
the Autism Genetic Resource Exchange (AGRE) Consortium 
(http://research.agre.org/program/descr.cffn). All AGRE-Pa- 
tients met diagnosis for 'autism' and each family was ascertained 
on the condition that at least two individuals were diagnosed with 
ASD [24], where diagnostic tests included autism diagnostic 
interview-revised (ADI-R) [25] and the autism diagnostic obser- 
vational schedule (ADOS) [26] . Regulatory review, approval, and 
oversight of AGRE's human subject research is provided by 
Western IRB (title: AUTISM GENETIC RESOURCE EX- 
CHANGE). In order to find new sequence variations with a 
putative function-altering effect, all exons and their flanking 
intronic sequences of CACNB2 were sequenced in 155 ASD 
patients. The number of sequenced cases enabled the detection of 
rare SNPs (f = 0.01) with a probability of 96% to detect it once and 
81 % to detect it twice. The power calculation was made according 
to Glatt et al. [27]. To assure a strong genetic background, 
patients from multiplex families with a severe phenotype were 
selected. 
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Figure 1. Splice scheme of human CACNB2 resulting in 9 splice variants of the Ca v |!2-subunit with the localization of the three 
mutations. Spliced exons are shown in light grey and conserved exons in dark grey. All nine splice variants express the mutation-carrying exon 5, 
while three of variants contain the localization of the third mutation in exon 7c [10]. 
doi:10.1371/journal.pone.0095579.g001 
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Table 1. Clinical Features of ASD patients with rare mutations in CACNB2. 





Description 


Patient 1 


Patient 2 


Patient 3 




(AU062603) 


(AU1 338304) 


(AU066204) 


Mutations 


c.334G>A 


C.42SOT 


c.553T>C 


(# NM_20 1597.2) 


p.C167S 


p.S197F 


p.F240L 


Genotype 


+- 


+- 


+- 


Gender 


M 


M 


F 


Epilepsy 


Yes 


Unknown 


No 


Type 


Complex partial 






Autism 


ADI-R 


Autism 


Autism 


Autism 


ADOS 


Autism 


Autism 


Autism 


Brain Imaging 


MRI 


normal 


Unknown 


Normal 



ADI-R, Autism Diagnostic Interview-Revised; ADOS, Diagnostic Observation Schedule-Generic; MRI, magnetic resonance imaging. 
doi:1 0.1 371 /journal.pone.0095579.t001 



The control group consisting of 259 matching controls was 
kindly provided by Dr. Heusch and Dr. Brodde, published in [28] . 
Based on the gender-dependent frequency of ASD [1], we chose 
an appropriate male to female ratio for both groups of about 5:1. 

Genotypic and DNA Sequence Analyses 

Oligonucleotides to all known exons of the CACNB2 gene were 
designed according to genomic sequences found in the Celera data 
base using Primer3. PGR amplification of DNA samples was 
carried out with GoTaq (Promega) according to the manufactur- 
er's protocol; annealing temperature was chosen according to 
Primer3. Mutational analyses were carried out with Mutation 
Surveyor. PGR fragments were purified using Exol and SAP 
(NEB), and sequencing was performed with an ABI 3730 
automated DNA sequencer. 

DNA Constructs and Site-directed Mutagenesis 

For functional analysis, mutations p.G167S and p.S197F were 
introduced in human p 2d cDNA (NM_201596.2) and the mutation 
p.F240L was introduced in P 2dE7c (NM_20 1597.2) by site- 
directed mutagenesis (Stratagene QuikChange Kit) and verified 
by sequencing. 

EGFP was used as reporter gene, which was coexpressed together 
with the P2-subunit by the bicistronic pIRES2-EGFP vector 
(Clontech). G167S forward primer 5'-gatagggcgattggtaaaagaaagctgt- 
gaaatcggattc-3'; G167S reverse primer 5'-gaatccgatttcacagctttcttt- 
taccaatcgccctatc-3'; S197F forward primer 5'-cagagagccaagcaa- 
gggaaattctacttcagtaaatcaggaggaaattcatcatcc-3'; S197F reverse 
primer 5 ' -ggatgatgaatttcctcctgatttactgaagtagaatttcccttgcttggctctc- 
tg-3'; F240L forward primer 5'-ggaaaactgcaggcttgctttggcggtttac- 
taca; F240L reverse primer 5'-tgtagtaaaccgccaaagcaagcctg- 
cagttttcc. 

Cell Culture 

In brief, HEK-293otlc cells stably expressing Ca v 1.2 subunit 
(GeneBank accession $JMM_000719) cloned from human heart 
[29] were grown in Petri dishes in Dulbeccos modified Eagle 
medium (PAA, Germany) supplemented with 10% FBS (PAA, 
Germany), penicillin (10 U/ml), and streptomycin (10 ug/ml; 
PAA, Germany). Cells were selected by geneticin (G418) (PAA, 
Germany) at a final concentration of 500 ug/ml. Cells were 



routinely passaged twice a week and incubated at 37°C under 6% 
C0 2 growth conditions. HEK-293 otic cells were transfected using 
SuperFect reagent (Qiagen). For whole-cell recordings, cells were 
transfected with a 2:3 ratio of either human WT or mutant P 2d - 
subunit and human a 2 5 1 -subunit [30]. 

Electrophysiology 

Whole-cell recordings in EGFP-positive cells were obtained 48- 
72 h after transfection. Immediately prior to recording, cells kept 
in 35-mm culture dishes were washed at room temperature (19- 
23°C) with bath solution. For constructs based on human 
P 2d cDNA (NM_201596.2) the bath solution contained (in mM) 
20 BaCl 2 , 1 MgCl 2 , 10 HEPES, 40 TEA-C1, 10 dextrose, and 65 
CsCl (pH 7.2 with TEA-OH) and pipette solution (in mM) 105 
CsCl, 25 TEA-C1, 11 EGTA, and 10 HEPES (pH 7.2 with 
CsOH). Holding potential was — 100 mV. For constructs based on 
human P2d_E7c cDNA (NM_20 1597.2) the bath solution 
contained (in mM) 5 BaCl 2 , 1 MgCl 2 , 10 HEPES, 40 TEA-C1, 
10 dextrose, and 92 CsCl (pH 7.2 with TEA-OH). Patch pipettes 
made from borosilicate glass (1.7 mm diameter and 0.283 mm 
wall thickness, Hilgenberg GmbH, Malsfeld, Germany) were 
pulled using a Sutter Instrument P-97 horizontal puller and fire- 
polished using a Narishige MF-83 microforge (Narishige Scientific 
Instrument Lab, Tokyo, Japan). Pipette resistance was 2-4 Mfl. 

Currents were elicited by applying test potentials of —40 mV to 
+60 mV from a holding potential of — 1 00 mV using Clampex 
software pClamp 5.5 and an Axopatch ID amplifier (Axon 
Instruments, Foster City, CA, USA). Voltage dependence of Ba~ 
current inactivation in HEK cells was determined with a two-pulse 
protocol. For P 2d cDNA (NM_20 1596.2) constructs the condition- 
ing first pulse (stepped from — 1 00 mV to +40mV) was held for 5s 
and the second pulse (of +10mV) for 125 ms; for p 2d _ E7(: cDNA 
(NM_20 1597.2) constructs the conditioning pulse was stepped 
from —75 mV to +5 mV for 1 s and the second pulse (of + 10 mV) 
was held for 125 ms. For both steady-state inactivation protocols 
the relative magnitude of inward current elicited during the second 
pulse was plotted as a function of the voltage of the conditioning 
first pulse. 

Data were analyzed using pCLAMP6 (Axon Instruments) and 
GraphPad 5 Prism. Currents were filtered at 2 kHz. Data were 
fitted to a Boltzmann function to obtain the half point (V0.5) and 
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(A) Patient 1(AU062603) Patient 2(AU 1338304) Patient 3(AU066204) 

C.334G>A; p.G167S C.4250T; p.S197F c.553T>C; p.F240L 
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Figure 2. (a) Chromatogram showing point mutations for patient 1-3. (b) The pedigrees show the occurrence of the disease phenotype and the 
inheritance of missense mutations in CACNB2. Circles and squares indicate females and males, respectively. Filled symbols denote diagnosis of autism; 
partially filled symbols denote individuals that are diagnosed within a broad spectrum of pervasive developmental disorders, (c) Alignments show 
conservation among different species for p.G167S, p.S197F and p.F240L suggesting an importance for the topology and function. 
doi:1 0.1 371 /journal.pone.0095579.g002 



slope factor (dV) for the voltage dependence of inactivation; Fits 
were performed after subtracting the offset from the peak values of 
the steady-state inactivation data. Offset was defined as a deviation 
from zero at the end of a fully inactivating current (after a 5000ms 
prepulse followed by a 125ms test pulse at +10mV). For voltage 
dependence of activation data were fitted by combined Ohm and 
Boltzmann relation I(V) = {V — Vr) x — it- according 

(l + cxp^^J 

to [31]. 

Statistical Analysis 

Student's unpaired two-tailed t-test (GraphPad Prism 5, 
GraphPad Software, San Diego, CA) was used to compare 
electrophysiological parameters gained with the several CACNB2 
mutations and their corresponding wildtypes. P<0.05 was 
considered significant. All p-values are listed in table S2. Data 
are presented as mean ± SEM. 



Results 

We performed a screening for mutations in CAC.NB2 by 
sequencing all 19 exons and their flanking intronic regions from 
155 patients with ASD. All variants found are listed in table SI. 
We filtered for variants causing non-synonymous substitutions, 
frameshifts, or splice site changes and identified six missense 
mutations. Since two of them occurred at similar frequencies in 
our patient cohort, our control cohort, and the NCBI database of 
short genetic variations (dbSNP), they were not further considered. 
One substitution ("Newl 1" in tabl SI.) is likely to be benign 
according to the prediction of polymorphism phenotyping tool 
(PolyPhen2); also the amino acid exchange of the minor allele is 
prevalent in other species (data not shown). Three non-synony- 
mous mutations were located in conserved domains of the protein 
(Fig. 1) and subsequently analyzed for functional effects. They 
were detected in three different patients who all were heterozygous 
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Figure 3. 1-V relationships. Currents were elicited from -40 to +60 mV in 10 mV increments with 5 mM or 20 mM Ba 2+ , respectively. The holding 
potential was - 1 00 mV. Current density of the variants p.GI 67S (N = 1 1 ) and p.SI 97F (N = 8) in (a) and p.F240L (N = 5) in (b) were compared with 
their respective WTs (|32d_WT: N = 9 and (32dE7c_WT: N = 5). 
doi:1 0.1 371 /journal.pone.0095579.g003 



carriers (Fig. 2a). The mutations were absent from 259 healthy 
controls. 

Patient 1 showed a variation at cDNA position 334 from 
guanine to adenine resulting in an amino acid exchange from 
glycine to serine (p.G167S, referring to NM_20 1596.2). The 
substitution is not annotated in dbSNP and according to the 
PolyPhen2 prediction 'probably damaging' for the CayP2 
function. Patient 2 had the cytosine at position 425 changed to 
thymine resulting in a substitution of phenylalanine for serine 
(p.S197F). This mutation is an annotated rare variant and shows a 
'possibly damaging' PolyPhen2 score. The third affected patient 
had the transition C.T553C which results in a phenylalanine-to- 
leucine change (p.F240L, referring to NM_201597.2). The 
missense mutation is not annotated and PolyPhen2 predicts a 
'benign' phenotyping score. The mutations p.G167S and p.S197F 
were found by analyzing the fifth exon, which is expressed in all 
splice isoforms; p.F240L was found in the alternative exon 7c 
(nomenclature according to Foell et al. [32]). All patients were 
diagnosed with autism by both ADI-R and ADOS. The clinical 
features of the patients are compiled in table 1 . 

The subsequently performed cosegregation analysis showed a 
low penetrance. Two variants showed incomplete segregation 
since some siblings, who were not obviously affected, were carriers 
of the mutation (Fig. 2b). Since complex-genetic disorders manifest 
depending on multiple factors, a low penetrance and incomplete 
cosegregation can be expected. In the literature, it has been 
suggested that pathogenic variants might be masked by epigenetic 
marks and thus are less penetrant [33] or that pathogenic variants 
follow an oligogenic pattern of inheritance, resulting in more 
severe or different phenotypes or are modified by the presence of 
additional variants [34]. Accordingly, unaffected carriers of a 
mutation might be possibly influenced by unknown genetic and 
epigenetic mechanisms [35] or might be subclinically affected. 

As members of the membrane-associated guanylate kinase 
(MAGUK) protein family [8-10], (3-subunits are comprised of two 
highly conserved central (GK and SH3) domains, flanked and 
interspersed by more variable N- and C-termini and a central 
HOOK domain (Fig. 1). The [S2-subunit is subject to extensive 
alternative splicing of the N-terminus and within the central 
HOOK domain, comprising the alternative exons 7a-c (Fig. 1) 
[36]. All three mutations are located within the HOOK domain 



(Fig. 1), which is known to modulate calcium channel inactivation 
[37,38]. The wildtype sequence of all three missense mutations 
revealed full conservation across species (Fig. 2c) indicating 
functional relevance. 

To investigate the effect of the mutations on calcium channel 
function, we heterologously expressed WT and mutant forms of 
the human [32 together with human a25l and EGFP in HEK-293 
cells stably expressing a human Cavl-2 pore-forming subunit, 
followed by whole-cell patch-clamp recordings. There are nine 
splice-variants of Ca v P2, with different inactivation kinetics [4]. 
The mutations p.G167S and p.S197F are present in all these 
Ca v fS2 splice forms, and there is no predominant splice form in 
the brain. In accordance with our hypothesis and deductions from 
the TS-channel model, we expect a reduction of the channels 
inactivation rate. To evaluate this effect in the mutations we 
performed whole-cell experiments using Ba + as charge carrier, 
since the TS-mutation selectively decelerates time-dependent 
inactivation and spares or even accelerates calcium-dependent 
inactivation [39]. For the quantification of mutants' influence on 
channel inactivation behavior we chose the [S2d splice variant as a 
suitable reference isoform, because of its relatively fast inactivation 
kinetics [4,8]. We compared the WT with the mutations p.G167S 
and p.S197F and the mutation p.F240L was analyzed separately 
with an appropriate WTE7c including its alternative exon 7c. The 
respective splice forms |32d (NM_20 1596.2) and (52dE7c 
(NM_20 1597.2) are - according to the GNF Expression Adas 1 
and 2 - expressed in fetal brain and in various adult brain regions. 

Peak current densities tended to be enhanced in the mutants 
P2d_G167S and |32d_S197F compared to P2d_WT. These 
mutations also caused a slight, non-significant shift of the I-V 
relationships towards more positive test potentials (Fig. 3a). 
Conversely, for the mutation of the P2dE7c_WT 
(P2dE7c_F240L), the I-V relationship revealed a trend towards 
enhanced channel activity at negative test potentials (Fig. 3b). 

Steady-state inactivation was examined using two-pulse proto- 
cols. Data and Boltzmann fits are depicted in Fig. 4a and 4b. The 
voltage of half-maximal inactivation (V0.5i nact ) did not differ 
significandy between the three mutant channel complexes and 
their respective controls (Fig. 4c). However, a significant flattening 
of the Boltzmann curves was observed for P2d_G167S and 
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Figure 4. Voltage-dependent steady-state inactivation (a, b) of Ba currents through L-type calcium channels. The ASD mutants 
p.G167S (N = 4) and p.S197F (N = 6) showed a significantly flattened slope of voltage-dependent inactivation compared to (32d_WT (N = 7). The third 
mutation (32dE7c_F240L. (N = 5) did not obviously differ from its corresponding (32dE7c_WT (N = 2). Half-inactivation potentials (V0.5 inact ) (c) and the 
slope factors dV (d) were obtained from the fits of individual experiments using the Boltzmann equation and averaging the results. Asterisk (*) marks 
a statistical significance (P<0.05) compared to the respective WT. 
doi:10.1371/journal.pone.0095579.g004 



(32d_S197F, as indicated by the slope factor dV (Fig. 4d). Numeric 
values are compiled in table 2. 

As described previously for the Timothy Syndrome, the most 
notable differences between mutants and WT were observed in the 
time course of whole-cell Ba 2+ currents. The mutant |32d_G167S 
and P2d_S1967F subunit induced a decelerated time-dependent 
inactivation of the channel complex within a 150ms test pulse 
(Fig. 5a). The WT with the alternative exon 7c inactivated more 
slowly than the WT with exon 7a. It was therefore examined using 
an extended test pulse duration (see scale bars in Fig. 5a). 
Compared with the adequate P2dE7c_WT the mutant 
P2dE7c_F240L showed an accelerated time-dependent inactiva- 
tion. The mutants and the WT with Exon 7a were recorded with 
20mM Ba i+ (compare [9]) and those with exon 7c were recorded 
with 5mM Ba _+ because of insufficient voltage-control with 



analyzed the extent of time-dependent inactivation as the 
percentage of current that has inactivated after 150ms of 
depolarization (% inactivation). As expected from the representa- 
tive whole-cell current traces, the currents inactivated significandy 
more slowly for all test potentials in case of P2d_S197F and at — 
20mV in case of P2d_G167S (Fig. 5b). In contrast, the mutant 
P2d_F240L led to faster, more pronounced inactivation compared 
to its corresponding WT (Fig. 5c). After 1000ms the currents of 
P2d_WT, P2d_G167S and P2d_S197F showed a comparable 
extent of inactivation of 94 ± 1%, 99 ±2% and 96 ±3% at OmV test 
potential. 



Discussion 

The complex genetic principle of origin underlying ASD 



i still 
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20 mM Ba 2+ and prolonged test pulse protocols (500ms). We to be elucidated [40] . So far, genome-wide association and linkage 
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Table 2. Current and voltage parameters (mean ± SEM) of the constructs tested. 







-P2d 


p2d_WT 


p2d_G167S 


p2d_S197F 


p2dE7c_WT 


p2dE7c_F240L 




N = 3 


N = 9 


N = 11 


N = 8 


N = 5 


N = 5 


Peak Current Density [pA/pF] 


-1.2 


-56.5 


-72.3 


-68.8 


-17.8 


-21.1 




±1.0 


±8.6 


±10.3 


±9.7 


±5.8 


±3.6 


V0.5 act [mV] 




-11.9 


-9.8 


-6.3 


-13.6 


-18.7 






±2.0 


±1.2 


±2.5 


±0.36 


±2.5 


Steady-state Inactivation 




N = 7 


N=4 


N=e 


N = 2 


N = 5 


V0.5 inact [mV] 




-38.7 


-47.9 


-35.7 


-33.3 


-28.2 






±1.8 


±5.6 


±3.7 


±0.7 


±1.9 


dV 




9.7 


16.9 * 


13.5 * 


8.1 


8.8 






±1.2 


±2.9 


±1.0 


±1.5 


±0.7 



V0.5 act , Voltage of half-maximal activation; V0.5 inac „ half-maximal inactivation voltage, dV, slope factor; -p2d, mock transfected cells; 
indicates P<0.05 versus corresponding WT. 
doi:1 0.1 371 /joumal.pone.0095579.t002 



studies presented inconsistent loci, reflecting a broad etiological 
heterogeneity and suggesting the influence of rare variants 
weighted by common susceptibility alleles. To put the spectrum 
and the genes into a pathophysiological context, an oligogenic 
model with epistasis has been assumed [41,42]. It describes - based 
on their level of biological function - a combination of multiple 
interacting genes resulting in ASD phenotypes. 

Until now, few potential pathophysiological mechanisms have 
been postulated [43,44]. The role of impaired calcium channel 
inactivation for ASD was revealed by Splawski and colleagues with 
the TS-mutation (CACNA1C: p.G406R) [2]. In more general 
terms, calcium signaling was found by a gene pathway analysis 
based on the 'Kyoto Encyclopedia of Genes and Genomes' 
(KEGG) as one of the three most significant pathways for ASD 
[45]. The concept is further supported by the association of 
mutations in various calcium channel genes with non-idiopathic 
ASD [2,13-16]. The activity and regulation of all L-type and some 
non-L-type calcium channels are dependent on their auxiliary (3- 
subunit; therefore, the highly differentiated P-subunit function is 
expected to be crucial for proper function of the calcium channel 
complex and homeostasis in the nervous system. 

Here, we present for the first time rare substitutions in highly 
conserved residues of the calcium channel subunit f>2 in three out 
of 155 individuals with ASD, all of which were absent from 259 
matching controls. According to the exome variant server (http:/ / 
evs.gs.washington.edu/EVS), 100 variations (35 of which are 
missense, splice or frame-shift mutations) are known for European 
Americans in CACMCB2. Thus, a statistical proof of an elevated 
frequency of rare variants in the CAGNB2 gene in ASD patients 
was considered not feasible. We did not test for differences in the 
mutation frequency between ASD patients and controls. Instead, 
we searched for most promising candidates of putative causal 
variants. Subsequently, by using patch clamp, the variants were 
tested for putative impairment of the function of the channel 
complex. Electrophysiological analyses in HEK-cells by whole-cell 
patch-clamp recordings demonstrate that all three missense 
mutations significantly alter the kinetics of the currents carried 
by the Ca v l .2 complex. Two mutations in the Ca v P2-subunit lead 
to deceleration of time-dependent inactivation of Ba 2+ currents as 
well as altered sensitivity of voltage-dependent inactivation. The 
third mutation shows a non-significant hyperpolarizing shift in 
current-voltage relation and an accelerated time-dependent 
inactivation. Of note decelerated time-dependent inactivation 



and incomplete voltage-dependent inactivation behavior are the 
biophysical hallmarks of the TS-mutation p.G406R in CACNA1C 
[2,17,39]. 

We observed a low penetrance of the three mutations in the 
ASD families under investigation, indicating the influence of other 
factors on the full expression of the condition. However, rare 
variants might contribute to the complex genetics and clinical 
heterogeneity of ASD. For comparison, an analogous sequencing 
approach studying CACMA1H as a candidate gene for ASD in 461 
patients revealed six non-synonymous mutations in conserved 
domains, all showing a low penetrance and incomplete segregation 
[15]. Expression of complex genetic disorders depends on multiple 
factors, therefore a low penetrance of causal variants is quite 
common. Unaffected carriers might be subclinically affected, other 
risk factors might have contributed to the phenotype in the 
affected individuals, or the described mutations may only act as 
modifiers of the phenotype. Even for the Timothy syndrome - one 
of the most penetrant monogenic forms of autism - ASD has a 
penetrance of ~75% only [46]. 

The putative effects of the Ca v P2-mutations are ample as the P- 
subunit also plays an essential role as an interaction partner with 
intracellular signaling machineries. For instance, Krey et al. have 
recently demonstrated a new pathway involving the TS-mutation 
which leads to a calcium channel activity-dependent and 
conformation-dependent dendrite retraction via Gem-induced 
RhoA-signaling [18]. Here, the Ca v P-subunits might be involved 
in two ways: firstiy they are known to regulate the channel activity, 
but secondly the interaction of Cay 1.2 with Gem, mediated 
through the P-subunit is important for the ability of the channel to 
control activity-dependent dendritic arborization [18]. In summa- 
ry, the role of the CavP2-subunit as a risk factor for ASD can not 
only be attributed to its interaction with diverse ASD-associated 
pore-forming Ca v Q(l-subunits, but also because the Ca v P2- 
subunit acts as a signaling hub and can link together different 
ASD pathways. 

The mutations presented here appear to follow a similar but 
milder mechanism of action that occurs in TS. TS presents with a 
multi-organ dysfunction, possibly indicating the consequences of a 
dramatic elevation of the intracellular calcium concentration. 
Mutations resulting in an even larger effect would likely abort the 
organism's development [47]. Thus the mutation of the TS can be 
viewed as an extreme of the viable spectrum of mutations within 
the calcium signaling pathway. Compared to this, the mutations 
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Figure 5. Time-dependent inactivation. Representative whole-cell Ba + current traces (a) of Ca v 1 .2/ot281 co-expressed with f32d_WT -subunit or 
co-expressed with the ASD mutants (32d_G167S, )32d_S197F or the (32dE7c_WT or its variant (32dE7c_F240L The currents were elicited from a holding 
potential of -100 mV by 150 ms (or 500ms) step depolarization of -40 to +50 mV with 20 mM Ba 2+ (for (32d variants) or 5 mM Ba 2+ (for Fi2dE7c 
variants), respectively. Analysis of the extent of time-dependent inactivation. (b, c) % Inactivation was analyzed as the remaining fraction of whole- 
cell current that has not inactivated after 150 ms of depolarization. |32d_WT (N = 12), (32d_G167S (N = 8), (32d_S1 97F (N = 1 1 ) (c) Similar analysis for 
P2dE7c variants p2dE7c_WT (N = 5), (32dE7c_F240L (N = 6). Asterisk (*) marks a statistical significance (P<0.05) compared to the respective WT, (**) 
marks P<0.01. 

doi:1 0.1 371 /journal.pone.0095579.g005 



P2d_G167S, P2d_S197F and P2d_F240L exhibited rather mod- 
erate effects on channel gating that nonetheless might suffice to 
unbalance neuronal calcium channel function. 

Because private and rare mutations seem to play a role in the 
predisposition to ASD [44], the discovery of rare variants with 
putative functional relevance might contribute to our understand- 
ing of the disorder's etiology. Integrating the data from TS [2] , the 
meta-analysis of psychiatric disorders [22], and functional studies 
on auxiliary P-subunits, we propose that inappropriate function of 
different components of the voltage-gated calcium channel 



complex can result in or may contribute to autism spectrum 
disorder. More detailed biophysical and cell-biological studies 
under physiological conditions are warranted for all such 
mutations. 

Supporting Information 

Table SI All found variants in the exons and flanking 
intronic regions of CACNB2 in autistic patients. 

(DOC) 



PLOS ONE | www.plosone.org 



8 



April 2014 | Volume 9 | Issue 4 | e95579 



CaVp2-Subunit Mutations in ASD 



Table S2 Current and voltage parameters (mean ± 
SEM) of the constructs tested. P-values stem from 
individual Student t-tests of the indicated channels. 

(DOC) 
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